Exact Density-Functionals for Correlated Dynamics on a Quantum Ring 
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We explicitly construct the time-dependent Kohn-Sham potential as a functional of the density 
and the initial state for a quantum ring with two interacting particles in a singlet state. In this 
case the initial state is completely determined by the initial density, the initial time-derivative of 
the density and a single integer that characterizes the (angular) momentum of the system. We 
give an exact expression for the exchange-correlation potential that produces the exact density of 
a non-interacting reference system with a different initial state. This is used to demonstrate that 
the Kohn-Sham procedure correctly predicts the density of a reference system without the need of 
solving the reference system's Schrodinger equation. We further construct the exchange-correlation 
potential and frequency-dependent exchange-correlation kernel for an analytically solvable system 
of two electrons on a quantum ring with a squared cosine two-body interaction. 

PACS numbers: 31.15.ee, 71. 15. Mb, 31.10,+z 



Time-dependent density- functional theory (TDDFT) 
[l], Q allows for an exact description of a many-body 
system in terms of an effective non-interacting system, 
known as the Kohn-Sham (KS) system. The effective 
external potential (known as the KS potential) in the 
non-interacting system is a functional of the density in 
such a way that the KS system has exactly the same 
density as the reference system. The essential compo- 
nent of the KS potential is the exchange-correlation (xc) 
potential that contains all nontrivial many-body effects. 
It depends on the initial states of the interacting and 
the KS system (initial-state dependence) as well as the 
density at all previous times (memory). Both features 
of the xc potential are, however, not well understood 
and consequently virtually all commonly used approxi- 
mations neglect them, which in important cases (dou- 
bly excited states, molecular dissociation, charge transfer 
etc.) can lead to large errors in the calculated properties 
[H, [2] • It is therefore highly desirable to have exact an- 
alytical functionals available for model systems that can 
serve as benchmarks and which can provide insight into 
the way memory and initial state dependence can be in- 
corporated into approximate functionals for real systems. 
In this work we explicitly construct such exact function- 
als of the density and the initial state for the case of a 
quantum ring (QR) with two interacting particles in a 
singlet state. We further give an analytic expression for 
the xc potential that produces the exact density of a non- 
interacting reference system with a different initial state. 
We use this to demonstrate that the KS procedure cor- 
rectly predicts the density of a reference system without 
the need of solving the Schrodinger equation of the corre- 
sponding reference system. We further consider the an- 
alytically solvable system of two electrons on a QR with 
a squared cosine two-particle interaction, for which we 
construct the xc potential and the frequency-dependent 
xc kernel of linear response TDDFT. 
The Runge-Gross theorem 3] and its generalizations 
d, § tell us that, for systems with general two-body in- 



teractions, the wave function is a functional of the one- 
particle density n(r,t) and of the initial state l^o), i- c - 
|*(t)) = |*([*o, n],t)). This also applies to the external 
potential that produces this wave function by propaga- 
tion of the initial state, i.e. v(r,t) = v([i>o,n],r,t). Thus 
the knowledge of the density n(r, t) is sufficient, in prin- 
ciple, to calculate all observables of the quantum system. 
A numerically efficient way to calculate the density is 
the KS construction, where one uses an auxiliary non- 
interacting quantum system (the KS system) that pro- 
duces the same density as the corresponding reference 
system of interest. The existence of such a KS system 
in the time-dependent case has been shown in [6| and 
more general proofs have since been put forward [J, [5| . 
Since the non- interacting system with initial state |$o) 
also obeys the Runge-Gross theorem, there is a unique 
potential v s (r,t) = v s ([$o, n], r, t) for this system that 
produces the density n. We can now define the xc poten- 
tial by the expression 



Here the Hartree potential is defined as 

v H ([n],r,t) = { dr'n(r',t)w(r,r'), 



(1) 



where w(r,r') is the two-particle interaction. We further 
define the KS potential by 

v KS [^ ,^ ,n\ = Wext + v n [n] + v xc [V ,<$>o,n] (2) 

where w cx t is the external potential in a reference sys- 
tem that, in general, is taken to be interacting. Sup- 
pose we use vks as the potential in a non-interacting 
system with initial state |$o) an( i we solve the resulting 
nonlinear Schrodinger equation (i.e. the KS equations) 
self-consistently for the density. Then for the density n 
that we obtain we must have, by definition of v s , that 
«Ks[*Oi^o,^] = V s [$o 5 «]- This then implies, by Eqs.([T]) 
and that w[5'o,n] = u ex t, and therefore the density 
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that we obtain is identical to the density obtained by 
solving the Schrodinger equation of the interacting sys- 
tem with initial state |^o) an( A external potential v cxt - 

We will now explicitly construct the functional 
v s [<I>o,n] and the functional WxcI^O) n] for the case 
that the reference system is non-interacting. This will 
be done for a one-dimensional two-electron system with 
periodic boundary conditions, also known as a QR. We 
start by making the Ansatz 



®(x,y,t) = <p(x,t)<p(y,t) 



(3) 



for the spatial part of the singlet state KS wave function. 
The full space-spin dependence is obtained by multiply- 
ing by the usual anti-symmetric singlet spin-function. 
We start by constructing ip as a functional of the den- 
sity. This orbital is determined by the one-dimensional 
Schrodinger equation 



id t (p(x,t) = 



2 x 



v s (x,t) (f(x,t), 



(4) 



with periodic boundary conditions on the interval [0, L] 
and starting from the initial state ipo(x) — ip(x,to). We 
then rewrite the orbital in terms of real- valued functions 
\tp\ and S as [J 0,1 

<p(x,t) = \tp(x,t)\exp(iS(x,t)). 

The periodic boundary conditions on the orbital then 
correspond to periodic boundary conditions on \ip\ and 
for the function S we find 

S(L, t) = S(0,t)+2irm (5) 
d x S(L.t) = d x S(0,t), (6) 

for some integer to. We see that the initial orbital 
ipo{%) = \fo( x )\ exp(iSo(x)) determines the choice of to 
since Sq(x) = S(x,to) must obey condition ([5]). The den- 
sity and current of the KS system are given by 

n(x,t) = 2\<p(x,t)\ 2 (7) 
j(x,t) = n(x,t)d x S(x,t). (8) 

which are connected via the continuity equation 

-d x j(x,t) = -d x [n(x,t)d x S(x,t)} = d t n(x,t), (9) 

and which expresses the local conservation of particles. 
This is a Sturm-Liouville equation Q parametrically de- 
pending on the time t and thus the density determines 
the phase function S(x,t) for a given set of bound- 
ary conditions (0 and ([6]). More precisely S is deter- 
mined uniquely up to a purely time-dependent constant 
since the constant function is eigenfunction of the Sturm- 
Liouville operator in © with eigenvalue zero and also 
satisfies the boundary conditions ([5]) and (j6]). Physically 
this freedom amounts to gauge freedom in the potential. 
Following similar derivations as in Q we then find that 



S([m, n], x, t) 



dyK t (x,y)d t n(y,t) (10) 
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where we defined 

K t ([n},x,y) = he{y-x)-e{x-y)] f 

2 J y n[z,t) 

T](x,t)r](y,t) 

fL dy ' 
JO n(y,t) 



with 6 the Heaviside function and 



r)(xt) 



dy 



n(y,t) 



dy 



n(y,t) 



All these functions are defined within the interval [0, L\ 
and can be extended periodically outside of it. At t = 
to this equation determines So(x) in terms of n(x,to), 
dtn{x,to) and m and therefore for a given choice of m 
the density completely determines the initial state tpo(x). 
Therefore, if we restrict ourselves to the product Ansatz 
([3]), there is for a given time-dependent density only a 
countably infinite number of initial states possible. 

We can deduce from Eqs. (0 and (fTOj) that the different 
values of m correspond to different total momenta P(t) = 
dxj(x, t) (or angular momenta if we interpret x as an 
angular variable ). Thus, the initial states which all have 
the same initial density and time-derivative of the density 
but different boundary conditions on the phase, describe 
differently rotating systems. 

The external potential v s can readily be expressed in 
terms of the orbital by inverting the Schrodinger Eq. ^ 
and we find 0, 0, 1 



v s ([cp],x,t) = 

l%\<p(x,t)\ 



id t (p{x,t) + \dlip.{x,t) 



<p{x, i) 



(11) 



2 ^(»,*)-i(^(*,t)) 9 

.ra^ + ^^ i 

\<p{x,t)\ \f{x,t)\ 2 



The last term on the right hand side vanishes as a con- 
sequence of the continuity Eq.® and we thus find using 
\<p\ = \J n/2 that 



JO n(z,t) 



(\ l 4\ 1 9 x y/n{x,i) 
v s ([m,n\,x,t) =- — - dtb([m,n\,x,t) 

z \/n{x,t) 

-±(d x S([m,n},x,t)) 2 , (12) 

which gives v s as a functional of n and the initial state 
(characterized by to). The potential v s [to, n] exists when- 
ever we have a unique S[m, n], i.e. for n > and 

the integrability conditions J Q dx \l/n(x, t)\ < oo and 

J Q L dx\d t n(x, t)\ < oo are fulfilled [5j. Thus we have an- 
alytically defined a density-potential mapping which is 
also explicitly initial-state dependent. We stress that the 
periodic boundary conditions on the wave function were 
essential in deriving Eq. (TT2")) . This excludes, for instance, 
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the example of a homogeneous electric field on a ring of 
constant density given in reference @ [l5| . 

Let us now consider the xc potential. The general def- 
inition of this potential is according to Eq. (fTJ) 

u xc [*o,$o,n] = v s [$ ,ri\ - vn[n] - u[*o,n]. 

The construction of this functional therefore requires the 
knowledge of the functional v[$>o,n}. This functional is 
not explicitly known. However, if the reference system 
is non-interacting then i^Oi ?i] = ^sf^O) ra ] an d since in 
that case v-n[n] = we find that 

Uxc[#o,$o,n] = v s [$ ,n] -u s [* ,n] (13) 

and hence the KS potential of Eq. ([2]) is given by 

v KS [® ,$> ,ri\ = Wcxt + v s [<S>o,n] - V s [9 ,n\. (14) 

Since for the case of our QR we know the functional 
w s [$0: n ] for the case of a product Ansatz initial state, 
we can apply expressions ffT3]) and (fT4")l . The xc poten- 
tial is then given by 

u xc ([m,m', n],x,t) = v s ([m',n],x,t) - v„([m,n],x,t) 

„ , „ / So7^t}\ . 2tt 2 (to 2 - m' 2 ) 
= 27r(m -m)d t \ rL K J + — 

\J0 n(z,t)J [n( X ,t) J Q ^yj 

2n(m-m')d x S{[0,n],x,t) (15) 
f L dz n ( x t) ' ^ > 

J0 n(z,t) V ' 

where d x S([0, n], x, t) is defined only in terms of n and 
d t n and corresponds to the spatial derivative of the first 
term on the right hand side of Eq. (flQ|) . For the ring 
system the integer m plays the role of the initial state 
^>o whereas m! plays the role of the state <f>o. The corre- 
sponding KS equations are thus 

id t ip(x,t) = ( - i<9 2 + v cxt (x,t) 

+ v X c([m,m',n],x,tfjcp(x,t), (16) 
n(x,t) = 2|(^(x,t)| 2 , (17) 

with tp(x,to) = tp™ (x). This equation now determines 
the density n(x,t) of the reference system when we pre- 
scribe Wext- We note that the xc potential is given 
only in terms of n and dtn. In contrast, the functional 
v s ([m,n],x,t) of Eq. (I12[) that reproduces a prescribed 
density via propagation of the KS equation also con- 
tains a second-order time-derivative of the density (in 
the term dtS as can be seen with the help of Eq. (1101) ). 
We therefore can explicitly see that the second order 
time-derivative of the density vanishes if we connect the 
two systems. This is an important fact which sometimes 
is overlooked in the literature and can lead to misun- 
derstandings about the KS approach [l], [l(| [H|. Thus 
the KS equations (flBl) and (fT?} have to be solved self- 
consistently. The initial state of the reference system 



tells what the initial density and initial time-derivative 
of the density is. The initial state of the KS system must 
then be chosen such that it gives these quantities cor- 
rectly. Since the KS potential in Ea. ([THl) only requires 
the knowledge of the density and its first time derivative, 
these initial data determine the solution uniquely. 

In the case that the reference system is interacting, we 
do not know v xc [^0: n] but we can still calculate it as 
a function of space and time for a specific density when 
we can solve the interacting system. This we will now 
do for the case of two electrons interacting via a squared 
cosine potential on a QR of length L. We start by the 
explicit construction of all the eigenstates of this system. 
The Hamiltonian is given (in atomic units) by 

H= -^+d 2 )+Acos 2 (|(x-y)), (18) 

where A is the strength of the interaction and x and y 
are the spatial coordinates for the particles along the 
ring. In this example v cxt is therefore simply the zero 
potential. Thus for repulsive interactions the particles 
will preferably be located where the interaction is zero, 
i.e. on opposite sides of the ring. For a two-electron 
system the eigenstates can be written as the product of 
a spatial wave function ^(x, y) and an (anti)-symmetric 
spin-function. We have a spin-singlet (spin-triplet) con- 
figuration if ^(x, y) is (anti)-symmetric with respect to 
an interchange of x and y, i.e 

*(x,y) = ±*( y> x) (19) 

where + refers to the singlet state and — to the triplet 
state. We further have the periodic boundary conditions 

*(x + L,y) = 

V{x,y + L) = *(x,y), 

with similar conditions on the spatial derivatives. It is 
convenient to introduce the center-of-mass R = (x + y)/2 
as well as the relative coordinate r = x — y. In terms of 
these coordinates the Hamiltonian of Eq. (fTS"]) attains the 
form 

H= -^-3 2 + Acos 2 (^). (20) 

The eigenstates Q(R, r) — ^(x,y) in the new coordinates 
then satisfy the equivalent property of (TT9l 

$(R,r) = ±$(JZ,-r) (21) 

and the periodic boundary conditions 

$(R + ^,r±L) = <i>(R,r) (22) 

and similarly for the spatial derivatives. With the Ansatz 
$(i?, r) — f(R)g(r) the Schrodinger equation can be sep- 
arated. The periodic boundary conditions on / and g 
become 

g(r + L) = ±g(r), (23) 
f(R + L/2) = ±f(R), (24) 
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and similarly for the spatial derivatives, where the signs 
on the right hand side of these equations must be the 
same for / and g in order to fulfill Eq. (|22|) . The equa- 
tion for the center-of-mass coordinate R becomes a free 
particle Schrodinger equation 

- \dlf{R) = ef(R) 

which has the eigenstates (up to normalization) 

f(R)=exp I — — I 

where the boundary conditions with ± in Eq. (|24p corre- 
spond to k being even and odd respectively. The energy 
eigenvalue is e = (kn/L) 2 . After changing coordinates 
to z — m/L the Schrodinger equation in the relative 
coordinate becomes 

d 2 z M{z) + [a-2q cos (2z)] M(z) = 0, (25) 

where we defined M(z) = g(Lz/ir). We further defined 

L 2 ( \\ 

XL 2 
q = 4**> 

with E the eigenenergy of the full Hamiltonian of Eq. (TT8l) . 
The boundary condition (1231) then becomes M(z + tt) = 
±M(zY Eq. ([25]) is the well-known Mathieu equa- 
tion [13]. The solutions are given by the Mathieu-sine 
and Mathieu-cosine functions denoted by SE(l, q, z) and 
CE(l, q, z) where / is a non-negative integer labeling cer- 
tain discrete values ai for the constant a in Eq. (l25l) . 
In the limit A — > (non-interacting case) we simply 
have CE(l, 0, z) = cos(Zz) and SE(l, 0, z) — sin(Zz) and 
ai = I 2 . We thus see that the ± signs in the bound- 
ary conditions (|2"5)) correspond to the case that I is even 
and odd respectively. From Eq. (f2"Tj) we see that the 
singlet and triplet case corresponds to the symmetry 
g(r) — ±g(—r) or equivalently M(z) = ±M(—z) for the 
Mathieu functions. This means that the singlet solu- 
tion corresponds to the Mathieu-cosine function and the 
triplet to the Mathieu-sine function. The full solution of 
the problem is therefore given by 

*+ (x, y) =N+ exp (j-k{x + y)^j CE (l, q,j(x- y)) 

*m(x,j/) =iV,-exp (^k(x + y)j SE (l,q,^-(x - y)) 

where + and — refer to the singlet and triplet cases re- 
spectively and 2Vi is a normalization factor. In both 
cases k and I need to be both even or both odd. The 
associated energy eigenvalues are 

E£ l =(iy[k 2 + af(q) + 2q], 



where a ; (q) are the characteristic values for the Mathieu- 
cosine and Mathieu-sine function respectively [13|. For 
q ^ the characteristic values obey (q) < aj" (q) < 
a t(q) < a 2 (?) ^ ■•■ 1 wnne in the non-interacting case 
af (0) = aj~(0) = I 2 . We thus nicely see how the two- 
particle interaction splits the degeneracy of the spin- 
singlet and spin-triplet states. Independent of the inter- 
action strength the ground state of the QR is the spin- 
singlet state ^oq(x, y). For large values of q the Mathieu 
functions become localized around z = tt/2 (and hence 
r = L/2) corresponding to the strongly correlated limit 
of well-localized electrons on opposite parts of the ring. 
The limit L — > 00 corresponds to q — > 00 and corresponds 
to a limit where the density goes to zero. This limit cor- 
responds to the famous Wigner crystal (T3 |. 
We can now construct the exact xc potential for a spe- 
cific density n that corresponds to a solution of the time- 
dependent Schrodinger equation with Hamiltonian (|18|) . 
For such a density the xc potential is given by 



v xc = v s [m',n] - va[n] 



since in our example ?; ex t = 0. The xc potential can 
be further split into an exchange (x) and a correlation 
(c) part i> xc = Vx + v c where, for our two-electron system, 
the x potential is simply given by v x = —(1/2)vh We 
choose the density n to come from a freely propagating 
superposition of two normalized eigenstates of our QR 



${x,y,t) = C y+ {x,y)exp(-iE+(t-t )) 
+ Ci*&(x,y) exp(-LE+(t - to)), 



which is a solution to the time-dependent Schrodinger 
equation. This wave function is properly normalized 
whenever Cg + Cf = 1. Note that both eigenstates have 
a constant density. If the constant Co is almost 1 (or 
0), the density of the system only deviates slightly from 
being homogeneous. If we look at a small QR, e.g. L = 1 
and different interaction strengths A, we find that even 
for small deviations from homogeneity the c potential is 
at least of the same order of magnitude as the x poten- 
tial. In this case, increasing the density variations by 
changing C'o make the correlation potential v c the domi- 
nant contribution to i> xc . A notable exception is an initial 
KS state that has approximately the right initial angular 
momentum (in the case of A = 100 and L = 1 this is the 
state ml = 1 as can be seen in Fig. (JTJ) ) . 
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b.) x potential 




o o 

c.) c potential (m'=0) 






d.) c potential (m'=1) 








FIG. 1: (color online). The density, x potential and c poten- 
tials for Co = 0.5 and m' — as well as m' = 1 (A = 100, 
L = 1). Note the change of scale between m' = and m' = 1. 

Here the c potential plus the x potential mainly needs 
to cancel the Hartree potential. The KS orbital would 
travel around the ring in approximately the right manner 
without any external perturbations. Besides the initial- 
state dependence one also clearly sees the non-locality of 
the c potential in time and space, as it has in general 
no obvious simple relation to the local density (see the c 
potential for to' = in Fig. (JTJ). If we go to larger QRs, 
e.g. L — 2n, the x potential becomes the dominant con- 
tribution to the xc potential. This seems counterintuitive 
since for this case the value of q is larger, corresponding 
to a more correlated state. It should, however, be remem- 
bered that the relation between the density profile (and 
hence the shape of v s ) and the electronic correlations is 
rather indirect. For example, the ground state density 
and KS potential of the QR are spatially constant, inde- 
pendent of the interaction strength. To get more insight 
into the influence of interactions, it is therefore more use- 
ful to study a two-point function. We therefore now will 
construct the xc kernel, which is defined to be the first 
functional derivative of u xc with respect to the density n, 
evaluated at the ground state density. 
The xc kernel is the central object of interest in linear- 
response TDDFT from which one can determine the per- 
turbative dynamics of the quantum system and its exci- 
tation energies. We start by calculating how the ground 
state spin-density reacts to small external perturbations, 
i.e. 

Sn(xa, w) — J dx'x(xa,x'a\oj)Sv(xa',oj) (26) 

a' 

(see e.g. in Refs. flU), where 



X(xa,x'cr',uj) = 



-(E P M -E )+ie 



J kl 

(*o|n(xV)|^)(^|n(xa)|* c 



with e > an infinitesimal, h(xa) the usual spin-density 
operator and — oo < k < oo and < I < oo (k and I are 
always either both even or both odd). Here with p = — 
we refer to the triplet state with spin function (<5 C r,t < ^T'4 + 
5<ji ^5 a ,_i)/V2 since the spin functions orthogonal to this 
one give a zero contribution in the sum. In a first step 
we can deduce using the periodicity of the solutions that 

(y \h(xa)\y+) = exp(i2Trkx/L)D + (k,l), 
(*o|n(a:o-)|^) = exp (i2nkx/ L) D~ (k,l){S u - 5 ia ), 

where 



D+{k,l) = N+N+ I drCE (o,q,^ CE (l,q,^-r 



x exp 



(-ijfcr) , 

D~ (k, I) = N+N t - / dr CE (o, q, ^r) SE (l, q, ^r) 

(-ii*) . 



x exp 



We note that the Mathieu cosine and sine are real and 
thus we have D ± (k, I)* = D ± (—k, I). After some manip- 
ulations of the general expression for the linear-response 
kernel we end up with 

x(xa,x'a',uj) = (u)5 aa , 

k 

+ /^( w )[i-<w]Kfc(z)a(z')*, (27) 



where 



Ck(x) 



exp (2s£«) 

7T~ 



2L (E±-E Q )\D±{k,l)\ 2 
r {u + ie) 2 -(E±-E ) 2 ' 

where the sum runs over all even values of I if k is even 
and over all odd values if k is odd. In the non-interacting 
case we find due to \D ± (k,l)\ 2 — > 5\ k ^i/(2L 2 ) the simple 
expression 



k 2 



^ {UJ) ~l{lJ ( w + ie ) 2 -i(-) 4 ,4 

and (u>) = 0. In linear-response (spin) TDDFT the 
interacting response function x is expressed in terms of 
the non- interacting one xvn i- e - 



X = Xo + Xo/hxcX, 

where the Hxc kernel is defined as 

p -l -l 

/Hxc = Xo _ X ! 



(28) 
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and summation as well as integration over reoccurring 
spin-space variables is implied. With the inverse kernels 
of Eq. (J2TJ) we find that 



t4M 



(w) 4vZ(u)v k (u) t 
4l/ fc ( W K M J 



The xc kernel is then trivially found by subtracting the 
interaction potential, i.e. / xc = /h xc — w. The restriction 
to k 7^ in the sum is a consequence of the fact that 
the response functions are only invertible in the space 
of functions orthogonal to the constant function, since 
a constant potential variation gives no density change. 
One can therefore add an arbitrary constant to the xc 
kernel in Eq. ([28|) without changing x- We have now fully 
characterized the behavior of the interacting particles on 
a QR in terms of the KS system for weak external per- 
turbations. The xc kernel exhibits a strong frequency de- 
pendence as it needs to shift the poles of \o and generate 
new poles in order to have the correct density response 



of the correlated system. Needless to say, this is strongly 
dependent on the interaction strength A and size of the 
QR L. 

In conclusion, we presented the exact form of the KS 
potential for a spin-singlet KS wave function and the ex- 
plicit form of the xc kernel of linear-response TDDFT in 
the case of a one-dimensional quantum ring with two in- 
teracting electrons. We have illustrated the capability of 
the presented results to serve as a new benchmark-system 
for approximations used in TDDFT and presented an ex- 
act KS construction. These results will be valuable in de- 
ducing new and more reliable functional approximations 
in TDDFT. 
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